Phases of the excitonic condensate in two-layer graphene 
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Two graphene monolayers that are oppositely charged and placed close to each other are consid- 
ered. Taking into account valley and spin degeneracy of electrons we analyze the symmetry of the 
excitonic insulator states in such a system and build a phase diagram that takes into account the 
effect of the symmetry breaking due to the external in-plane magnetic field and the carrier density 
imbalance between the layers. 
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I. INTRODUCTION. 

The excitonic insulator— ~— was predicted theoretically 
four decades ago in 3D semiconductors and then in spa- 
tially separated layers of electrons and holes.^ Since 
then, an excitonic insulator has been searched for in a 
variety of systems. The excitonic insulator is a mate- 
rial where the electron-hole excitonic correlations lead 
to the formation of a gapped state characterized by the 
order parameter resembling a superfluid condensate of 
excitons. Such a correlated state has been observed in 
double-quantum well semiconductor structures in quan- 
tizing magnetic fields£~— After the experimental discov- 
ery of graphene^— it has been discussed as a possible 
candidate for experimental realization of the excitonic 
insulator stated— sparking the on-going debated— 
about the critical temperature T c of the excitonic conden- 
sate transition in a two-layer graphene system. Various 
estimations of T c for such a system lay in a wide region 
of magnitudes from milli-Kelvms 2 ^' 27 up to Kelvins 2 ^— 
and further up to the room temperature! 22 ' 23 The consid- 
ered system22 r 23 i 25 ~ — consists of two parallel, separately 
controlled graphene monolayers, in which external gates 



(a) i 2 _ (b U 

, Zeeman splitting, 2e z 

GATE . , . GATE 



GRAPHENE 




FIG. 1. (a) The excitonic condensation due to an electron- 
hole pairing is studied in the system of two spatially separated 
graphene monolayers with an excess of electrons on layer 1 and 
a lack of electrons on layer 2. (b) The schematic phase dia- 
gram of the excitonic condensation in the system at different 
values of a Zeeman splitting and different values of the asym- 
metry between Fermi energies in layer 1 and 2, ez — M-B|h i s 
the Zeeman energy in an in-plane magnetic field h. 



induce a finite density of electrons in the layer 1 and holes 
in the layer 2, Fig. [IJa). Recently the two-layer graphene 
system has been obtained experimentally^ - — 

In this paper, we extend the existing theory of the ex- 
citonic insulator state in a two-layer graphene system: 
we analyze a symmetry of the excitonic insulator and 
classify its phases. As a result a phase diagram of the 
excitonic insulator is built, that takes into account the 
effect of the symmetry breaking due to the Zeeman split- 
ting and the asymmetry between electron/hole densities 
in the layer 1 and 2. A phase diagram, Fig. [T|b), contains 
3 phases: B,B' and Ai. Transitions between phases are 
found to be of the first order. These transitions are sub- 
ject to the use of an in-plane magnetic field and a varia- 
tion of external gate voltages, leading to different charge 
carriers densities in layers: the density of all electrons 
m e in layer 1 (which corresponds to the Fermi energy 

Ep = hvy/Tmi e /2), and density of holes n^h in layer 2 
(which corresponds to the negative Fermi energy in the 

-hv^/Trn 2 h/2). 



layer 2, Ef 



The B phase, Fig. [ljb), exists when there is no mag- 
netic field and charge carriers densities are the same in 
both layers n\ e — ri2h (ke. when Ep = \Ep |). The 
B' phase exists at the same condition n\ e = nib, but 
when an in-plane magnetic field is applied, which causes 
a Zeeman splitting of energies of electrons with different 
spin projections. The Ai phase exists when a symme- 
try of charge carriers density is violated, e.g. n\ e > n^hi 
and when the corresponding splitting of Fermi energies 

E 



(i) 



(2) 

Ep | is equal to the Zeeman splitting due to an 
in-plane magnetic field. 

The diversity of obtained phases is due to the high 
symmetry of the normal ground state, which can be bro- 
ken in several different ways leading to a variety of phases 
possessing different symmetry groups. A well known ex- 
ample of the system with the diversity of phases due to 
various normal state symmetry breaking is liquid Helium- 
3 1??"" 42 In liquid Helium-3 the symmetry of the order pa- 
rameter can be changed by correspondent external pa- 
rameters, leading to phase transitions. 

The analysis in this paper is organised as follows. Sec- 
tion |TT] describes the theoretical model of the considered 
system. Pairing of electrons and holes within mean field 
theory is introduced in the Section IIIII Section IIVI pro- 
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vides symmetry analysis and the phase classification of 

the excitonic correlated state. Section IVl contains de- 1>CjP) s 2, £ , p, — s 

tailed description of the most symmetric phases, their \ / 

properties are summarizes in Table |TJ Results are dis- V 
cussed in Section [VlJ 



II. TWO-LAYER HAMILTONIAN 



Graphene^r— is a gapless semiconductor with the 
Fermi surface consisting of two distinct points, K + and 
K_, called valleys. Near these Fermi points electrons 
have a linear dispersion E(p) — ±vp, with a velocity 
v R3 10 8 cm/sec^ here p — |p|, p = k K± is the mo- 
mentum of an electron relative to the Fermi point. Using 
external gates, one can independently tune the carrier 
density in each of the two graphene flakes.— Neglect- 
ing tunneling, the electrons in the two layer graphene 
system initially can be described with the Hamiltonian 
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-H\2, here the single particle 



part of the Hamiltonian is 



E 

Z,C,P,« 



{svp-Ef) a\ 



(1) 



the operators aj ^ s (ai.(. p . s ) create (annihilate) an elec- 
tron on the I = 1,2 layer on the ,s = +/— conduction 
or valence band with momentum p = p(cos </> p , sin p ), 



E { p> and E%> are the Fermi energies, which correspond 
to charge carrier densities in the layers. The index £ de- 
notes 4 different pairs of spin projection (f, j,) and valleys 
(K + , K_). In the Hamiltonian flayer the terms H\\ and 
H22 take into account the intra-layer interaction. These 
terms can be ignored in the following studies, provided 
that one uses a screened inter-layer interaction in the 
term H\2- Hence in H12 we keep only those terms that 
contribute to the BCS mean field theory^ absorbing 
other contributions into a renormalization of the velocity 
and the Fermi energy in the single particle part ([lj of the 
Hamiltonian 



,(2) 



H 12 = - £ nip-p , |) 1 + S5,C ° S(0P ~ V) 



p,p , S.S 

X E ! 'u p ~ 4,C',p',- S ' a U,PV' a 2,C',p,- s -( 2 ) 

C.C 

The scattering process, described by H12, is shown on 
Fig. [5J The function V(q) denotes a screened Coulomb 
interaction in the static limit V(q) — V(q,uj <C q). The 
factor [1 + ss' cos(4> p — </> P ')] /2 in Eq.([2|) reflects chiral 
properties of electrons related to the sublattice composi- 
tion of electronic Bloch wave functions! 44 ' 45 These chiral 
properties of electrons result in the suppressed backwards 
scattering if an electron does not change the energy band 
upon scattering (ss' = +), otherwise (ss' = — ) the elec- 
tron can not forward-scatter 42 



FIG. 2. A typical transition which is described by H12 in 
Eq. @. Indices I — 1 or 2, £, p, s denote a layer, a pair of the 
spin projection and valley, a momentum of an electron and a 
conduction s = + or valence s — — band. 



III. EXCITONIC PAIRING, MEAN FIELD 
ORDER PARAMETER 

The excitonic insulator state of the electron- hole liquid 
is characterized by the electron-hole correlations on the 
Fermi surface, Fig. 3. Mathematically it means that 
in the excitonic insulating state there is a non-vanishing 
ground state average F of electron operators 



Fccv(P) = (4,c',P,-s ai .C,P.+s) 



(3) 



For the existence of the non-zero ground state average 
F it is crucial that Fermi surfaces for electrons and holes 
coincide i 46 ' 48 Due to the electron-hole symmetry of the 
energy spectrum in graphene, the electron-hole excitonic 
correlations §3§ in the considered system are most devel- 
oped when the density of electrons in the layer 1 is equal 
to the density of holes in the layer 2, n\ e = U2h, or, in 
terms of Fermi energies E ( p ] = -E ( p\ FigEJ However 
apart from this condition there can be certain other ex- 
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E { p = -E^=v PF 
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FIG. 3. The excitonic electron-hole bound state in the two- 
layer graphene. The left hand side of the figure shows the 
electron's spectrum in graphene layer 1 and 2. An electron 
on the Fermi surface in the layer 1 is shown as a fulfilled 
circle. Absence of an electron on the Fermi surface in layer 2 
is shown as an empty circle. Closed line around both circles 
represents an excitonic pairing, which is developed due to a 
Coulomb interaction (shown as a wavy line). The right hand 
side of the figure shows the coincided Fermi circles in both 
layers at Fermi momentum pf- 
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ternal conditions when excitonic correlations ([3]) can be 
developed. Thus, although the excitonic insulator state 
disappears when the symmetry n\ e = n2h is violated by 
external gates, we show below that excitonic correlations 
can be restored by the in-plane magnetic field. Based on 
the detailed analysis of excitonic correlations in mono- 
layer graphene in the in-plane magnetic field, which is 
done by Aleiner and co-authors^, we show that the ex- 
citonic insulator state can exist in various phases in the 
two-layer graphene system. 

In order to study phases of the excitonic insulator state 
at different external conditions, firstly we apply the stan- 
dard mean-field approximation^ We assume that the 
product of the operators a\ >, _ s ai,f,p jS weakly deviates 
from its non-vanishing ground state average. We expand 
the interacting part -H12, Eq. ([2]), of the Hamiltonian 
flayer up to the linear order with respect to these small 
deviations and neglect constant terms. The mean field 
Hamiltonian of the system becomes 

H mi = H s . p . + ^2 K,Cp,s A CC':s(P) a 2,C',p,-s + H - c -L 

(4) 

where H.c. stands for "Hermitian conjugate", and 

A C c, s (p) = -E F cc^(p')nip-p'i) 

p',s' 

1 + SS' CQS{(/) P - (ftp/) 

x 2 ' (5j 

Quantities A^< jS (p) form the matrix A of the order pa- 
rameter. Index £ denotes 4 different pairs of spin pro- 
jections and valleys (f K+,t K-,4- K + ,| K_). Thus in 
the spinCgwalley space the order parameter is given by the 
4x4 matrix A with matrix elements given by (|SJ). For 
brevity we omit index s and momentum p in the notation 
for the order parameter A. 

For further analysis it is convenient to 
rewrite the Hamiltonian (j4| as follows: H m { = 
E Cp ,s*CP-- ff mf(p,s)* c , p , s , where # c , p , s 

(oi,C,p,+s) a2£,p,-s) T ) and 

tt f \ ( (svp-E F )l A \ 

H **(P>°) = [ A t -(svp-E F )l)' (6) 

Here all elements of the matrix H m { are 4x4 matrices 
in the spintgivalley space: diagonal elements have struc- 
ture of the identity matrix 1 in this space, whereas A is 
given by some 4x4 matrix, whose structure is identified 
in this paper for each phase of the excitonic correlated 
state. The matrix of the order parameter A describes 
the correlations between conduction/ valence electrons in 
the layers 1 and 2 below a critical temperature T c . 

Nevertheless, the phase classification can be made re- 
gardless of the value of the transition temperature T c . 



Assuming that the excitonic insulator state can be ob- 
served in two-layer graphene system, we analyze the sym- 
metry of the mean field Hamiltonian Q and the order 
parameter A. As a result the classification of all phases 
of the excitonic insulating state of the two-layer graphene 
system is presented in the next Section, and a detailed 
discussion of each phase is presented in Section |V| 



IV. SYMMETRY ANALYSIS OF THE 
CORRELATED STATE 

The analysis in this section is based on the idea of 
breaking of the initial symmetry of the hamiltonian by 
the order parameter. The initial symmetry group G 
of the Hamiltonian i^layer is formed by global unitary 
transformations of an electronic single-particle state in 
the 4-component spintgivalley space independently in the 
layer 1 and 2. These transformations are represented by 
independent matrices U^ 1 - 1 and U^ 2 - 1 in layer 1 and 2 re- 
spectively. Therefore the group G is given by the direct 
product of corresponding unitary groups U4 

G = ui 1] xUi 2) . (7) 

Unitary group uf\ I = 1,2, consists of 4 x 4 unitary 
matrices U") which perform transformations of electron's 
operators in the I— th layer as follows: 

ai,(,p,s->^2uf ] c a LC , p , s . (8) 

C 

Thus, as it is seen from Eqs. @ and ([5]), under symmetry 
transformations ([8]) the order parameter A transforms as: 

A^U (1)t AU (2) . (9) 

This implies that the Hamiltonian of the system is not 
invariant under the action of the group G any longer. 
However for any fixed non-zero A there is always some 
subgroup H of the group G, H C G, H ^ G, such that 
all transformations from the group H do not transform 
A, i.e. the order parameter A remains invariant: 

U^AU^A. (10) 

Such transformations U^ in layer 1 and U^ in layer 2 
form a symmetry group H 

( U fug>) effcG - 

Only transformations from the group H leave the ground 
state of the excitonic insulator invariant, i.e. only these 
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transformations leave the mean field Hamiltonian (j4]) and 
(|6"|) invariant: 

( o ug)tJ^CP,*)^ J u g ) J=^ f (p ! 4 

Thus the symmetry group G of the initial uncorrelated 
normal ground state of the system is broken down to the 
symmetry group H of the ground state of the excitonic 
insulator. 

All transformations from G, which are not included 
in H, form the factor-space G/H. These transforma- 
tions change the order parameter A, however they do 
not change the energy of the corresponding ground state. 
Therefore the manifold of all matrices A, which can be 
obtained by transformations from G/H, form a degener- 
acy space of the order parameter. Consequently the man- 
ifold of the correspondent ground states form a phase of 
the correlated state. It is important to notice, that all 
these ground states within the same phase are described 
by the same symmetry group H, which is a symmetry 
group of the order parameter. Therefore phases of a cor- 
related state can be classified by the symmetry group H 
and the degeneracy space of the order parameter. 

The phase classification presented in this paper is also 
reminiscent of the classification of the various degen- 
eracy spaces of the order parameter in liquid Helium- 
3,38 -40,42 This classification principle was used in the 
determination of superconducting phases in nontriv- 
ial superconductors^ and superfluid phases in liquid 
Helium-341^ 

In order to classify phases of the excitonic insulating 
state we are going to classify the degeneracy spaces of 
order parameters with the same symmetry. For this we 
consider the condition (fT0|) and use the method of a sin- 
gular value decomposition^ It allows us to represent an 
arbitrary matrix A as a product of a unitary matrix V* , 
a diagonal matrix D with real non-negative numbers on 
the diagonal, and another unitary matrix V. Applying 
the singular value decomposition to the order parameter 
at any given values of s and p we obtain 

A = Vt(p)D s (p)V s (p). (13) 

However, first of all, we notice that in the considered 
system, the lowest ground state energy is realized when 
matrices V, V do not depend on momentum p and index 
s = ±. The reason for this is that in such a case in the 
expression for the ground state energy there is a cancel- 
lation of the product of matrices V and (V and V^) 
into a unit matrix. Such a cancellation leads to the max- 
imal negative contribution to the ground state energy, 
and therefore the energy of the ground state achieves its 
minimal value. If we assume that unitary matrices in 
Eq. (|13p depend on momentum p and index s = ±, then 
unitary matrices at different momenta p, p' and differ- 
ent indices s, s' do not cancel each other, which increases 



the ground state energy comparatively to the previous 
case. Thus, we conclude, that in order to realize the low- 
est energy of the ground state, matrices V and V can 
not depend on the momentum p and index s. Therefore 
the singular value decomposition of the matrix of order 
parameter becomes 

A = V+DV, (14) 

where in the right hand side of the equation (I14[) only 
matrix D depends on p and s, but for brevity we omit 
these indices. 

Second, all transformations from the group G, includ- 
ing those from the factor-space G/H, do not change the 
diagonal elements of the matrix D, but change matrices 
V, V into any other unitary matrices. Thus if we intro- 
duce the notations 

V't = Tj(i)tvt, V' = VU< 2 \ (15) 

then under the transformation ([H])-® the order parame- 
ter will transform in the following way: 

A = V f DV — > A' = V' f DV'. (16) 

Thus under this transformation the diagonal matrix D 
does not change. Recall that the degeneracy space of the 
order parameter is obtained by acting on the order pa- 
rameter A by all transformations from the group G (here 
transformations from subgroup H will not change the 
order parameter while remaining transformations from 
factor-space G/H will create the degeneracy space of the 
order parameter). As long as only matrices V and V are 
changed by transformations from G, we obtain that the 
degeneracy space of the order parameter and a phase of 
the correlated state are determined only by the diagonal 
elements of the matrix D. 

Finally, from the condition (|TU| we have found that 
all possible degeneracy spaces of the order parameter are 
classified by numbers of equal and different diagonal el- 
ements in matrix D. In the case of physically relevant 
phases there are additional restrictions on the diagonal 
elements of matrix D. Thus among all possible matri- 
ces A only physically relevant order parameters satisfy 
the self-consistency equation. For the phase classifica- 
tion it is sufficient to consider the BCS self-consistency 
equation for the order parameter. Diagonalizing the self- 
consistency equation by unitary matrices from Eq. (|14l) . 
one obtains 4 equations for diagonal elements of the ma- 
trix D, each equation corresponds to some value of index 
£. These equations have the same structure and depend 
on the Fermi momentum pp. If the Fermi momentum pp 
is the same for all types of electrons (for all indices £), 
then these 4 self-consistency equations are identical, and 
apart from a trivial zero solution they have the same non- 
zero solution. Hence in such situation in physically rele- 
vant phases the arbitrary diagonal element in matrix D 
can be equal either to other non-zero diagonal elements, 
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TABLE I. The phase classification of the excitonic insulating state with respect to following two external parameters: (1) 
Zeeman energy ez = /is|h| in an in-plane magnetic field h, and (2) an asymmetry between the electron density ni e in the layer 
1 and the hole density n 2 h in the layer 2. 



External 
conditions 


Symmetry group G 
of the two — layer 
Hamiltonian 


Phase 


Matrix structure 
of the order 
parameter A 

(here V, V G U±) 


Symmetry group H 

oi tne 
order parameter A 


aim[G/-o J 


Single 
particle 
spectrum 


nie — n,2h 
ez = 


x U™ 


B 


V 




16 


gapped 


A' 


V 1 " Diag[l, 1,1,0] V 


x t/i 1 ' 2 ' x U[ 2) 


21 


gapless 


A[ 


V 1 " Diag[l, 1,0,0] V 


x t/ 2 (1 ' 2) x t/ 2 (2) 


20 


gapless 


A' 2 


V 1 " Diag[l, 0,0,0] V 


0™ x [/f ^ x 


13 


gapless 


nie — ri2h 
ez > 


t/ 2 (lt) x U[ u) x U { i f) x t/f w 


B' 


\ 


V v y 

',v G ?7 


2 




8 


gapped 


nie > n 2h ^ 
ez > 


Ui n) x x x tf 


Ai 




v G U 2 




[/< u) x t/ 2 (lt ' 2t) x t/ 2 (2i) 


4 


gapless 



a Here we assume parameters to be tuned so that a Fermi surface of electrons only with spin up in layer 1 coincides with the Fermi 
surface of holes with spin down in layer 2, for details see Fig. \E\ 



or be equal to a zero. The application of the in-plane 
magnetic field in principle changes such a description be- 
cause of Zeeman splitting. However, in the case where 
the Fermi energy is much greater than the Zeeman en- 
ergy, Ep ez, the magnetic field does not change the 
situation essentially as long as it is possible to neglect the 
difference between Fermi momenta for electrons with op- 
posite spin projections in the self-consistency equations. 
Therefore four self-consistency equations on four diagonal 
elements of the matrix D become approximately iden- 
tical also when relatively small in-plane magnetic field 
(ez <C Ep) is applied. 

The non-zero solution of these equations is given by 
the gap function g s (p) which at the Fermi surface (s = 
+, |p| = pf) determines a gap in the single-particle exci- 
tation spectrum. Thus we conclude that in all physically 
relevant phases the matrix D in the singular value decom- 
position (fl4|) of the order parameter A consists of zeros 
or non-negative diagonal elements which approximately 
are equal to the gap function g s {p)- 

Substituting the obtained result into Eq. (fl4|) . we ex- 
tract the gap function as a multiplier. Thus we con- 



clude that the order parameter A in all physically rele- 
vant phases has a form 

A^ 5s (p)V t DV. (17) 

Here the matrix D is a diagonal matrix with or 1 on the 
diagonal. The representation (|17[) becomes approximate 
in the case of the applied in-plane magnetic field with 
the condition ez <C Ep. The dependence of the order 
parameter A on variables s and p is completely given by 
the function g s {p)- 

Having the matrix of the order parameter provided, the 
symmetry group H is found from the equation (I10[) . For 
this the matrix of the order parameter is represented as a 
single value decomposition (fl"4|) . The constant matrices 

V and V are absorbed into matrices U^' and of 

the global symmetry transformations from the symmetry 

group H. Then the equation (|10[) connects two unitary 

(i) + ~ ( 2 )~ + 
matrices VU^ V T and W H V T and the diagonal matrix 

D with or 1 on the diagonal. Thus the matrices U# 

(2) 

and of transformations from the symmetry group H 



are obtained. 



V. PHASES 

In this section we provide a detailed description of 
phases of excitonic insulator state in two-layer graphene 
system. Results of this section are summarized in Table 
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A. The B phase. 

Firstly we consider the situation when there is no ex- 
ternal magnetic field and when the charge carrier densi- 
ties in layers are the same n\ e — n^h- In such case the 
symmetry group G of the two-layer Hamiltonian of the 
system in normal state is given in Eq. 0. Under the 
mentioned conditions the Fermi circle in the conduction 
band in layer 1 coincides with the Fermi circle in the va- 
lence band in layer 2 due to the electron-hole symmetry 
in graphene. Hence the non- vanishing ground state aver- 
age F, Eq. dU), can be formed by all species of electrons. 
Taking into account that the ground state with the lower 
energy is more stable, we consider the phase when ex- 
citonic correlations are developed among all species of 
electrons. In such cases the order parameter matrix A 
and the matrix D, Eq. (TIT)) , are not degenerate matrices. 
Moreover, because there is only one Fermi circle for all 
species of electrons, the most stable ground state is char- 
acterized by the matrix D in Eq. (1171) with equal non- 
zero diagonal elements, i.e. D is an identity matrix. As 
discussed in the previous section, this conclusion follows 
from the consideration of the self-consistency equation 
on the order parameter. Thus substituting D = 1 into 
the equation ([T7]) we obtain the following structure for 
the order parameter in spin<X>valley space: 



A = <? s (p)V, VeU 4 . 



(18) 



Such a structure of the order parameter determines the 
symmetry group H of the ground state and the degen- 
eracy space of the order parameter, and consequently it 
determines the phase of the excitonic insulator. 

The symmetry group H of the ground state in the con- 
sidered phase can be found as a group of all unitary trans- 
formations U# in layers 1 and 2, which leave the 



order parameter invariant, Eq. (1101) . Solving the condi- 
tion (fTQ)) with the order parameter (fT8|) we obtain ma- 



trices and of symmetry transformations in the 
layer 1 and 2 correspondingly. Thus, having TTjj and 



T (2) 



U 



(2) 
H ' 



ug>=u, 



u 



(2) 
H 



v f uv, 



(19) 



we can express an arbitrary element of the group H, Eq. 
(|11[) , which transforms electron operators in layers 1 and 



FIG. 4. In the B phase within the excitonic paired state the 
electron on the Fermi surface in the layer 1 (grey circle) is 
characterized by the index £ in the spin® valley basis $, and 
the absent electron (white circle) in the layer 2 is characterized 
by the same index £ but in the spin® valley basis V$, which 
is transformed by the matrix of the order parameter V. 



2 according to Eq. 



Omitting indices £, p, s we have: 



a i 

a-2 



U 

o vtuv 



a i 

a-2 



(20) 



Here the matrix V is given by the fixed matrix of the 
order parameter (|18[) . The unitary matrix U is present 
in transformations in both layers. This means that the 
symmetry group H of the ground state in the considered 
phase (|18p consists of the combined transformations in 
both layers 1 and 2. The unitary group of the combined 
transformations in layers 1 and 2 is denoted as Jj\ 



(1,2) 



U 

vtuv 



(1,2) _ 



= H. 



(21) 



r(i,2) 



The combined transformations from the group U\ 
can also be described in terms of generators of these 
transformations. For this, each element of the group is 
written as the exponential function of the element of the 
group's algebra 



U 

vtuv 



exp 



H 



e J7, 



(22) 



where 9 is a vector of real variables and a vector Th 
consists of generators of the group H. For the considered 
phase these generators are: 



(IV), 



A m 
VtA m V 



to = 0,1,.. .,15. (23) 



In contrast to Eq. (|2"3")) transformations, which change the 
order parameter and create a degeneracy space G/H, are 
described by the following generators: 



G/H) 



A m 
-VtA m V 



m = 0, 1, 



15. 



(24) 
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Here the 2x2 block matrices (r#) m and (Tq/ h ^ act 
in the space of layers 1 and 2, the matrix A m acts on the 
spin®valley basis in the layer 1 and matrices ±VtA m V 
acts on the basis $ in layer 2. The spin® valley basis 
$ is the same in both layers. Matrices X m are 4x4 
Hermitian traceless matrices of generators of transfor- 
mations from the unitary group U4,. The total number of 
generators Tq/jj, Eq. (jM]) . equals to the dimension of 
the degeneracy space G/H, dim [G/H]. In the B phase 
dim[G/H] = dim[G] - dim[H] = 32 - 16 = 16. 

The electron operators in the second layer can be trans- 
formed by the matrix of the order parameter: Va 2 — > a' 2 , 
see Eqs. (ITS)) and (gj, or, equivalently, the spin<x>valley 
basis in layer 2 can be transformed by the matrix of the 
order parameter. In such a case from Eqs. (|20|) and 
(|23f it follows, that the transformation from the group 
H, in contrast to the transformation from G/H, can be 
represented by identical transformations in both layers. 
These identical transformations act by the same matrix 
U on the spin<8>valley basis $ in the layer 1 and on the 
transformed spin® valley basis V<E> in the layer 2, Fig. 
21 Hence the matrix of the order parameter {T5J defines 
the relative unitary rotation of the spintgivalley basis $ 
in layer 2 with respect to layer 1. It signifies the relative 
symmetry breaking: the ground state is not invariant un- 
der unitary transformations of the basis $ in one layer 
relatively to the basis $ in another layer. The basis $ 
in layer 1 is "locked" relatively to the basis $ in layer 
2 by the matrix of the order parameter which defines 
the relative unitary rotation of one basis with respect to 
another. 

Because of the presence of the relative symmetry 
breaking by the order parameter (fT5)) the phase dis- 
cussed here resembles the superfluid B phase in the liquid 
Helium-3^-ii>^ 

In B phase the matrix of the order parameter is not 
degenerate. It means that all species of charge carriers 
develop excitonic correlations, therefore a single particle 
excitation spectrum is gapped. 

The external conditions (ez — 0, n\ e = n-ih) for the 
B phase can be violated by an in-plane magnetic field 
or by external gates. However the excitonic correlations 
continue to exist in the B phase until the difference of 
radiuses of Fermi circles is bigger than 2g + {pp)/v, where 
g+(j>p) is a gap in a single-particle excitation spectrum 
in the B phase. Indeed such behavior can be seen, if 
one creates an asymmetry between charge carriers den- 
sities in layers, which can be expressed in terms of a 
shift SEp > of Fermi enetgies: Ep^ = Ep + SEp, 

E F = —Ep + 8Ep. Substituting these values to 
the mean field Hamiltonian (j4]) and finding its eigen- 
values, one obtains^ two branches of excitation spec- 
trum £« (p) = \J {svp — Ep) 2 + g 2 {p) ± SEp. At values 
SEp = g+(pp) one of branches of excitation spectrum 
becomes zero at s = +,p = pp. At this situation the ex- 
citonic pairing stops being energetically favorable and the 
system appears in the normal state via a first order phase 



transition. In similar way, when Fermi circles for charge 
carriers with opposite spin projection are separated by 
2g+(pp)/v due to an in-plane magnetic field, excitonic 
correlations between charge carriers on these Fermi cir- 
cles vanish. This fact is schematically shown in the phase 
diagram FigJTjD: at the borders of the B phase in the 
phase diagram excitonic correlations are no longer ener- 
getically stable and the excitonic insulator state trans- 
forms into either a normal state or into another phase 
via a first order phase transition. 

B. The AqjA'^A^ phases. 

In this section we consider phases under the same ex- 
ternal conditions as in B phase, thus the symmetry group 
G is again given by Eq.Q. We consider phases where or- 
der parameters arc characterized by the degenerate ma- 
trices of rank r < 4. In such phases only a part of electron 
species develop excitonic correlations, therefore the single 
particle excitation spectrum is gapless for certain species 
of electrons. The matrix of the order parameter can be 
chosen as follows, compare with Eq.(fl"7|): 

A = . 9s (p)VtDiagM,c,0]V, \\VeU A . (25) 

Here the diagonal matrix Diag determines the order pa- 
rameter in phases, which are denoted as A' , A' x , A' 2 : 
numbers (a, b, c) are given by (1, 1, 1) in A' phase, (1,1,0) 
in A' x phase and (1,0,0) in A' 2 phase. Using the trans- 
formed electron operators Vai in the layer 1 and Va 2 in 
the layer 2, see Eqs. ((4]) and (|25]l . the self-consistency 
equation on the order parameter becomes diagonal, and 
only first r out of four equations for diagonal elements 
will have non-zero solutions. We assume that in the self- 
consistency equations we can use the screened interaction 
among charge carriers in the system in normal stated In 
this case self-consistency equations in A' Q , A! 1 , A' 2 phases 
are identical to self-consistency equations in B phase, 
therefore their non-zero solutions are given by the same 
gap function g s (p). 

Substituting the order parameter in the symmetry con- 
dition (|10p one obtains matrices of symmetry transforma- 
tions in layer 1 and 2, for example for A\ phase one gets 

<^t(» u °,)v, ug^Vt^^V, (26) 
where 

uedf' 2 ', u' ec/ 2 (1) , u"GC/ 2 (2) . (27) 

Here, similarly to the B phase, the 2x2 unitary matrix 
u determines the combined unitary rotations of the first 
two components of the spingivalley basis V$ in layer 1 
and the first two components of the spin£g)valley basis 
V<I> in layer 2. Therefore such phase is characterized by 
a partial relative symmetry breaking. Remaining matri- 
ces u',u" G U2 determine independent unitary rotations 
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of the other 2 components of corresponding spin(g)valley 
basis in layers. These other 2 components correspond to 
quasiparticle's states which do not contribute to the exci- 
tonic condensation, their single particle excitation spec- 
trum is gappless. Therefore only 2 out of 4 electron's 
species are involved in the excitonic condensation. Be- 
cause of this, such a phase is similar to the superfluid 
A\ phase of liquid Helium-3, where paired states with 
only one spin projection S z = +1 are present in the 
condensate^ The A\ phase of Helium-3 exists only in 
magnetic field, in order to underline the stability of the 
phase in the absence of the magnetic field we denote the 
phase discussed here by an additional prime, therefore 
it is denoted as the A[ phase of the excitonic insulator. 
Other phases with degenerate matrices of order parame- 
ters are denoted as A' and A' 2 . 

In phases A' , A[, A' 2 the number of non-zero diagonal 
elements in the diagonal matrix Diag determines the rank 
r of the symmetry group of combined unitary rotations, 

(1 2) 

denoted as Ur ■ Zeros in the diagonal of the matrix 
Diag correspond to electron states which do not develop 
excitonic correlations, and, therefore, these states can be 
unitary transformed independently in each layer. 

Consequently the symmetry group H for A' , A[, A' 2 
phases can be easily identified. For example, the symme- 
try group H for A\ phase is following: 



H = Ui 1] x U. 



(1,2) 



x U, 



(2) 



(28) 



The dimension of the degeneracy space is calculated as 
follows: for A' phase dim[G/H] = 32 - 1 - 9 - 1 = 21; 
for A[ phase dim[G/#] = 32 - 3 x 4 = 20; for A' 2 phase 
dim[G/H] = 32 - 9 - 1 - 9 = 13. 



B' phase 



(J) Electron with spin up 

(^Absent electron with spin down 
VPF kp v 





Pf - ez/v 




PF + ez/v 





FIG. 5. Excitonic correlations in the two-layer graphene 
system with an in-plane magnetic field h in the case of equal 
charge densities in layers rii e = ri2h- Because of a Zeeman 
splitting 2ez there are two Fermi circles with radiuses pf ± 
ez/v. The absent electron with a particular spin projection 
is considered as a quasiparticle (hole) with an opposite spin 
projection. 



a Zeeman energy ez should be bigger than a gap in the 
excitation spectrum in the B phase, tz > 9+{pf)- In 
such a case, due to the initial equality of charge carrier 
densities n± e = n 2 h in the B phase, the two Fermi cir- 
cles in layer 1 coincide with two Fermi circles in layer 2. 
Thus it leads to the appearance of two different Fermi 
circles in the system, Fig. [5] Consequently, the electron- 
hole pairs, which appear on different Fermi circles, have 
different properties: thus such electron-hole pairs have 
different spin projection, +1 or —1, Fig. [5] Also due to 
slightly different Fermi momenta, electron-hole pairs on 
different Fermi circles are characterized by slightly differ- 
ent gap functions. Hence in the spintgivalley basis $, 



C. The B' phase. 

In this section we consider the two-layer graphene sys- 
tem in an in-plane magnetic field. Our analysis is based 
on the comprehensive study by Aleiner and co-authors 19 
of the spontaneous symmetry breaking in graphene sub- 
jected to an in-plane magnetic field. When an in-plane 
magnetic field is applied, the Fermi circles for quasipar- 
ticles with different spin projections become separated 
due to a Zeeman splitting. Such splitting changes the 
symmetry group G, Eq. (|7|), of the initial Hamiltonian 
-f^layer toward a direct product of 4 unitary groups U2, 

G^fxC/W'xt/f'xJ/f*. 

Each of these U2 groups transforms a valley space of elec- 
trons with corresponding spin projections in one layer, 

(it) 

e.g. U 2 transforms electrons with spin up in layer 1. 

The B' phase can be obtained from the B phase by the 
application of an in-plane magnetic field. Such magnetic 
field should be big enough to break the excitonic correla- 
tions in the B phase and to split Fermi circles. Therefore 



«D = (tK + ,tK_,|K + ,|K_), (29) 

the order parameter has the following structure (compare 
with Eq. (PH) 

A=f ° (30) 

Here matrices v and v are unitary 2x2 matrices, which, 
by analogy with the B phase, determine the relative uni- 
tary rotation of a valley space of electron states with one 
spin projection in layer 2 relatively to electron states with 
another spin projection in layer 1. Functions g' s (p), g'J (p) 
are gap functions, which differ from each other only be- 
cause of the presence of a Zeeman splitting. However 
when the Fermi energy is much bigger than Zeeman en- 
ergy, Ep ^> tz-, the difference between these functions is 
negligible and they are approximately equal to the gap 
function in the B phase, g' s (p) ~ g" a (p) ~ <?s(p)- Thus, 

a-*<p)(2;). 
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Using this approximation, the symmetry group H of the 
order parameter can be found from the condition (jlOp . 

As a result one obtains matrices and Ujj of the 
transformations (| 1 1|) from the group H in the layer 1 
and 2 respectively (matrices are written in the basis (|2"9")l 
in each layer): 

«s>-(;§). < u (T;„)' < 3 » 

The corresponding electron operators are transformed as 
follows: 

■ f-i Jf. Q \ \ 

Oi,t -> uai.f, a.2,i -> v T uva 2 ,i, u € C/ 2 ' % (32) 

— ► uai.|, a2,t — >• v^uva2,t, u G E/g ' . (33) 

The unitary 2x2 matrix u, Eq. ([32"T) . determines a sub- 
group of the group H, which consists of combined unitary 
rotation of valley space of electrons with spin up in layer 
1 and electrons with spin down in layer 2: u G ZT^ . 
The unitary 2x2 matrix u, Eq. (1331) . defines another cor- 

responding subgroup of the group H, u G C H. 

Hence in the phase considered here the symmetry group 
H of the order parameter is given by direct product of 
two subgroups 

H = U^' 2i) xU^' 2t) . (34) 

Using expressions for groups G and H in the B'— phase, 
we found that the degeneracy space G/H is 8- 
dimensional, dim[G/H] = 4x4-4-4 = 8. It is also 
defined by the structure of the order parameter (|3"01) . i.e. 
here the degeneracy space is determined as a space of all 
possible unitary 2x2 matrices v and v. Because of the 
non-degenerate matrix of the order parameter, the single 
particle excitation spectrum in this phase is gapped. 

Similarly to the B phase, the excitonic correlations in 
the B' phase cease to exist when the external conditions 
(ez > g+(pF), Tiie = ri2h) are perturbed, i.e. when Fermi 
circles in different layers are separated for the energy in- 
terval which is bigger than a double value of a gap in the 
single-particle excitation spectrum. Thus, in particular, 
in the schematic phase diagram Fig. O^b) at the border 
of the B'— phase (when a symmetry n\ e = n^h of charge 
carriers densities is violated) the ground state of the sys- 
tem transforms to an uncorrelated normal ground state 
via the first order phase transition. 

D. The Ai phase. 

In contrast to B and B' phases, where all species of 
charge carriers develop excitonic correlations, in this sub- 
section we discuss another possible realization of the ex- 
citonic insulator state in the two-layer graphene system. 
We show that the excitonic correlated state can exist in 



A\ phase 




Pf + $Pf + ezfo 



FIG. 6. Excitonic correlations in the A\ phase. Starting from 
a zero magnetic field the asymmetry between charge carrier 
densities n\ e > niu is created. In terms of Fermi energies it 
means E^ > \Ep\. A magnitude of an in-plane magnetic 
field is chosen such that a Zeeman energy ez satisfies the 
condition: Ef' — ez = \Ef'\ + ez, where (Ef' — ez)/v is 
the radius of the Fermi circle for electrons with spin up in 
layer 1, and {\E { p\ + e z ) /v is a radius of the Fermi circle for 
electrons with spin up in layer 2. Both of these Fermi circles 
are situated at the same Fermi momentum pf- Therefore 
two out of four Fermi circles coincide, leading to excitonic 
correlations between only half of electron's species. Using the 
expression for the Fermi energies Ef^^ = ±vpF + vSpF, 
Pf > SpF, where Spf = (y/nie — v /n 2h)v / ^/2 5 the condition 
on the Zeeman energy is the following: ez = vSpF- 



the presence of an in-plane magnetic field and a specially 
chosen asymmetry in charge carriers densities in layers. 

In order to achieve a necessary external conditions, 
firstly we consider the two-layer system without a mag- 
netic field and with equal charge carrier densities in lay- 
ers. Under such conditions the spectrum of electrons in 
both layers have only one Fermi circle at the Fermi mo- 
mentum pf- By changing the external gate voltages we 
create asymmetry between charge carrier densities in lay- 
ers: ni e > n2h- Thus the Fermi circle in the layer 1 is sit- 
uated at the momentum pp + SpF, and the Fermi circle in 
the layer 2 is situated at the momentum pf — 5pf, where 
SpF > 0, and pf is the Fermi momentum in the case 
nie = Ti2h (he. in B and B' phases). It is assumed that 
the separation between Fermi circles is big enough to pre- 
vent the development of excitonic correlations. Keeping 
the chosen values of densities, we switch-on an in-plane 
magnetic field with such a magnitude that the Zeeman 
energy ez is equal to the energy shift of each Fermi sur- 
face, tz = vdpFi Fig. [5J The presence of an in-plane 
magnetic field signifies that the symmetry group G in 
such case is the same as in B' phase. 

The external conditions mentioned above lead to the 
situation when only two out of four Fermi circles coin- 
cide: both the Fermi circle of electrons with spin up in 
the layer 1 (a Fermi circle with a smaller radius in the 
layer 1) and the Fermi circle of electrons with spin up 
in the layer 2 (a Fermi circle with a bigger radius in the 
layer 2) are situated at the same Fermi momenta pf- 
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Thus electron-hole pairs are formed on these two coin- 
cided Fermi circles. Notice, that a total spin projection 
of such an electron- hole pair is equal to a zero in contrast 
to electron-hole pairs with spin projections +1 or —1 in 
the B' phase. Fermi circles for electrons with spin down 
in both layers do not coincide with any other Fermi sur- 
faces. Therefore corresponding electrons and holes are in 
a normal state (i.e. they do not participate in excitonic 
correlations), their single-particle excitation spectrum is 
gapless. 

Thus for the phase discussed here only half of all elec- 
tron species in the system develop excitonic correlations. 
It is reflected in the order parameter, whose structure 
in the basis (|29| in both layers is given by the following 
expression: 



Using the expressions for initial symmetry group G and 
Eq. (|4"0"j) . we find that the degeneracy space G/H in this 
phase is 4 dimensional, dim[G/H] = 16 — 3 x 4 = 4.. 
It is determined by the manifold of all possible matrices 
v G IJ2 in the structure of the order parameter ([35]). 

The subgroups and are present in both 

groups G and H, the appearance of the correlated state 
does not change them. Therefore for the phase discussed 
here the initial symmetry is broken only partially. Ac- 
cording to the similarities with the superfluid A\ phase 
of liquid Helium-3 (i.e. that the phase described here ex- 
ists in magnetic field and has a partial relative symmetry 
breaking ^ 41 ' 46 ), the phase discussed in this subsection 
was denoted as A\ phase. 



A = fc(p)(;[). (35) 

Here the gap function g s (p) is the same as in other 
phases due to the same Fermi momentum pp in the self- 
consistency equation and due to the approximation of 
the interaction among charge carriers in all phases by the 
screened interaction among charge carriers in the system 
in normal state. Similarly to the phases discussed previ- 
ously, the unitary 2x2 matrix v in the order parameter 
(|35p determines a relative unitary rotation of the valley 
space of electrons with spin up in layer 2 relatively elec- 
trons with spin up in layer 1. 

Solving the condition (|10[) with the order parameter 
(|35p we found that transformations from the group H 
are represented in layer 1 and 2 by following matrices 
and respectively (both matrices are written in 
the basis ([29]) in each layer): 

<={i°)- u M v, "»°»)- tm > 

Here matrices u, u',u" are unitary 2x2 matrices. The 
matrix u performs a combined unitary transformation of 
a valley space of electrons with spin up in both layers, 
in addition the valley space of electrons in layer 2 are 
rotated by the order parameter (|35|) . compare with Eq. 
(|32| . The valley space of electrons with spin down is 
transformed by the unitary matrix u' in layer 1 and by 
the unitary matrix u" in layer 2 correspondingly: 

—t uai,-)-, d2,t — > v^uva2,f 5 u € , (37) 

01,4. -> u'a u . u'e[/ 2 (u) , (38) 

a 2 ,4 -> u"a 24 , u"eu¥ l) . (39) 

Thus the group H consists of the direct product of 3 
unitary groups, 

H = U { ^ xU^' n) xU? l) . (40) 



VI. RESULTS AND DISCUSSIONS 

In the present paper we consider a two-layer graphenc 
system where external gate voltage induces a finite den- 
sity of electrons in one layer and holes in another. As- 
suming that the transition temperature T c towards exci- 
tonic insulator is high enough so that it can be observed, 
we classify phases of such correlated state. In order to 
obtain different excitonic correlations and therefore dif- 
ferent phases we propose to use parallel to graphene lay- 
ers magnetic field and perpendicular to graphene layers 
electric field. 

Firstly we consider the Hamiltonian of the two-layer 
graphene system. We recognize that the ground state is 
characterized by a high symmetry group - the group of 
unitary rotations of spintgivalley space of electrons in each 
layer independently. Below a transition temperature T c 
this symmetry is reduced by a non-zero order parameter 
towards a symmetry of the excitonic insulating ground 
state, which consists of electron-hole pairs with electrons 
on one layer and holes on another. Following the BCS 
theory of superconductivity, we identify the condition for 
such electron-hole pairing, determine the order parame- 
ter and build a BCS-like mean-field theory of the exci- 
tonic insulator. Analyzing a symmetry breaking of the 
initial ground state by the order parameter, we consider a 
condition that mutually determines the order parameter 
and the corresponding symmetry group of the excitonic 
insulator ground state. Using a singular value decom- 
position of the matrix of the order parameter, for each 
phase of the excitonic insulator we obtain a correspond- 
ing symmetry group of the ground state, a structure of 
the order parameter and its degeneracy space. The re- 
sults of a phase classification of the excitonic insulator are 
shown in TableQJ the most energetically stable phases are 
shown in the phase diagram, Fig. [ljb), and on Figs. [U 

It is important to notice that the excitonic correlations 
in all phases discussed in this paper origin from the co- 
incided Fermi surfaces at approximately the same Fermi 
momentum pp (we use assumption Ep 3> ez, where the 
Fermi energy Ep is determined in the system without a 
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magnetic field and with equal densities of charge carriers 
in layers n le = nzh)- Thus assuming that the interaction 
among charge carriers is the same in all phases (i.e. that 
the effect of excitonic correlations on the screening of the 
interaction can be neglected^!) , we obtain that the en- 
ergy gap in the single particle excitation spectrum in all 
phases is determined by the same self-consistency equa- 
tion. Therefore the transition temperature T c estimated 
from the self-consistency equation 4 ^ should be the same 
for all phases. 

At a temperature lower than the transition tempera- 
ture T c , transitions between phases in the phase diagram 
are found to be of the first order. Phases of excitonic in- 
sulator have different properties: thus the electron-hole 
pairs in B' phase have total spin projection +1 or — 1, 
Fig. [5j whereas in the Ai-phase a total spin projection 
of an electron- hole pair is equal to zero, Fig. [6] 

According to the number of estimations of the criti- 
cal temperature in the considered system, the most op- 
timistic estimation gives values of T c close to a room 
temperature^ However this estimation^ 2 - does not take 
screening of the Coulomb interaction into account, ex- 
plaining it by the assumption of the first order phase 
transition in the system. Some other estimation o 25 ' 27 
point on the improbability of observation of excitonic 
condensation due to extremely low transition tempera- 
ture < lmK, (T c « 1Q- 7 E F ). According to Refs^i 2 ! 
the reason for low transition temperature lays in the 



effective screening of the Coulomb interaction by a big 
number N of species of electrons! 25 ' 51 In the considered 
system N — 8, which is given by product of 2 val- 
leys, 2 spin projections and 2 layers. Such large N in- 
creases screening and makes excitonic condensation not 
so effective, as in the monolayer graphene- 5 - 2 " and espe- 
cially in the monolayer graphene in a magnetic field ) 19 ' 53 
where T c can reach value up to 10 -4 Ep. However, re- 
cent investigations, based on a detailed treatment of the 
screened Coulomb interaction 2 ^ 30 ' 54 ' 55 and on a con- 
sideration of a multi-band pairin g 31 ' 54 or pairing with 
nonzero momentum 5 - suggest that the transition tem- 
perature T c can be sufficiently big for the experimental 
observation of the excitonic insulator in the considered 
system. Together with recent experimental realization 
of the two-layer graphene system 3 ^— it provides a hope 
that the phase diagram of the excitonic insulator, Fig. 
Hlb), under favorable conditions 5 - 7 - - — will be observed ex- 
perimentally. 
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